Scaling of Local Slopes, Conservation Laws and Anomalous Roughening in Surface Growth 
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We argue that symmetries and conservation laws greatly restrict the form of the terms entering the long wave- 
length description of growth models exhibiting anomalous roughening. This is exploited to show by dynamic 
renormalization group arguments that intrinsic anomalous roughening cannot occur in local growth models. 
However some conserved dynamics may display super-roughening if a given type of terms are present. 
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Recent theoretical and experimental studies of self-affine 
kinetic roughening have uncovered a rich variety of novel fea- 
tures 1 1 ] . In particular, the existence of anomalous roughening 
has received much attention. Anomalous roughening refers 
to the observation that local and global surface fluctuations 
may have distinctly different scaling exponents. This leads 
to the existence of an independent local roughness exponent 
ai oc that characterizes the local interface fluctuations and dif- 
fers from the global roughness exponent a. More precisely, 
global fluctuations are measured by the global interface width, 
which for a system of total lateral size L scales according to 
the Family- Vicsek ansatz 1 2] as 



W{L,t)=lPf{L/t l l'), 
where the scaling function f(u) behaves as 

if u < 1 
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The roughness exponent a and the dynamic exponent z char- 
acterize the universality class of the model under study. The 
ratio j3 — a/z is the time exponent. In contrast, local sur- 
face fluctuations are given by either the height-height correla- 
tion function, G(l, t) = ([h(x + l,t) — h(x, t)] 2 }, where the 
average is calculated over all x (overline) and noise (brack- 
ets), or the local width, w(l,t) = (([h(x,t) - (ft,)/] 2 }/} 1 / 2 , 
where (•••}/ denotes an average over x in windows of size 
I. For growth processes in which an anomalous roughening 
takes place these functions scale as w(l,t) ~ \/ G(l, t) = 
^/aG/* 1//z )> w i tn an anomalous scaling function |3j|4j] given 
by 



fA{u) 



u ai "" if u < 1 
const, if u> 1 
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instead of Eq.(|2j. The standard self-affine Family-Vicsek 
scaling |2] is then recovered when a — <Xi oc . 

This singular phenomenon was first noticed in numerical 
simulations of both continuous and discrete models of ideal 
molecular beam epitaxial growth 

Anomalous roughening has later on been reported to occur in 



growth models in the presence of disorder mill, electro- 
chemical deposition models 1 14], chemical-vapor deposition 
lHHl . etc. 

Remarkably, anomalous roughening has also been reported 
in many experimental studies including molecular-beam epi- 
taxy of si/sum) m, sputter-deposition growth of Pt on 

glass I17|. Cu electrodeposition 1 18], growth of Fe-Cr super- 
lattices il9ll . propagating fracture cracks l2oll . fluid-flow in 
porous media fclll . tumor growth l l22ll . etc. 

Nowadays it has become clear |3l l23ll that anomalous ki- 
netic roughening is related to a non-trivial dynamics of the 
average surface gradient (local slope), so that ((V/1) 2 ) ~ t 2K . 
Anomalous scaling occurs whenever k > 0, which leads 
to the existence of a local roughness scaling with exponent 
Q/oc — ct — zk [23]. Also, it has recently been shown that 
the existence of power-law scaling of the correlation func- 
tions (i.e. scale invariance) does not determine a unique dy- 
namic scaling form of the correlation functions JU On the 
one hand, there are super-rough processes, a > 1, for which 
ctioc — 1 always. On the other hand, there are intrinsically 
anomalous roughened surfaces, for which the local roughness 
ai oc < 1 is actually an independent exponent and a may take 
values larger or smaller than one depending on the universal- 
ity class (see |4j, 12411 and references therein). The existence 
of intrinsically anomalous roughened surfaces is a most in- 
triguing observation and leads to the still open question con- 
cerning the basic physical features (symmetries, form of the 
nonlinearities, conservation laws, non-locality, etc) required 
for anomalous roughening to occur in surface growth. 

In this Letter we show that local models of surface growth 
driven by random noise cannot exhibit intrinsic anomalous 
roughening. Nevertheless, some models with conserved dy- 
namics may display super-roughening if the required terms 
are present (see below). Our results are based upon a dynamic 
renormalization group analysis of stochastic growth models 
that include only local interaction terms (derivatives of the 
height field) and the deposition noise is Gaussian and un- 
correlated in space and time. This framework is very broad 
and includes many universality classes of growth, in particu- 
lar stochastic equations used to theoretically describe molec- 
ular beam epitaxy as well as other types of thin-film growth. 



2 



Our results imply that disorder and/or non-local effects are re- 
sponsible for intrinsic anomalous roughening in experiments 
in different systems EUl . We also argue 

that the anomalous scaling exponents found in numerical sim- 
ulations of discrete growth models (see for instance 01) 
must be effective and corresponding to a non-universal tran- 
sient regime due to the use of too small and too short time 
scales compared with the true asymptotic regime. 

Scaling of the local slopes in surface growth.- We are in- 
terested here in surface growth models with local coupling 
among degrees of freedom, i.e, growth equations that include 
only derivatives of the height. In the long wavelength limit, 
we consider surface growth in d + 1 dimensions described by 
the Langevin equation 

dth = Q(yh)+ri(x,t), (4) 

where h{x, t) is the height of the interface at substrate posi- 
tion x and time t. The function Q(Vh) defines a particular 
model and incorporates the relevant symmetries and conser- 
vation laws. In particular, invariance under translation along 
the growth and substrate directions as well as invariance in 
the election of the time origin rule out an explicit dependence 
of Q on h, x and t. Growth is driven by an external noise 
ry(x, t), which represents the influx of particles in deposition 
processes. The noise is usually considered to be Gaussian, 
uncorrelated in space and time, and either non-conserved, 
{r}{x,t)r){x',t')) = 2D s ( x - - *')> or conserved, 

(r) c (xL,t)r) c (-x',t')) = -2D V 2 c5(x - x')S{t - t'), depend- 
ing on the model under consideration. 

The local slope fluctuations scale as ({Vh) 2 ) ~ t 2K . When 
k > the local slope fluctuations become a relevant scale (in 
the growth direction) and local scaling behavior is expected to 
be anomalous, Eq.Q, instead of standard self-affine Family- 
Vicsek. One can exploit this simple observation 1 23] to obtain 
the local roughness exponent ai oc = a — zk by analyzing the 
corresponding time evolution equation for the local derivative 
T = Vh 

d t T = Vg(T)+ry c (x,i), (5) 

where the noise r\ c = Vr] is now conserved. Assuming peri- 
odic boundary conditions in Eq.@ implies Vh = and one 
then finds that the global width of the interface T(x, t) is 

W r (t) = <[T(x,i) - Y^W) 1/2 = <(W> 1/2 , (6) 

which immediately leads to Wy{t) ~ t K 12311 . Therefore, 
by studying the scaling behavior of surface T(x, t) one can 
obtain the anomalous time growth exponent k and therefore 
the local roughness exponent of the surface h{x,t) through 
the scaling relation ai oc = a — zk. 

In all respects we can regard the slopes T(x, t) as a grow- 
ing surface on its own and study its roughening properties, 
if any. Specifically, time scale invariance of l|6) implies 
that the global width is expected to obey Eqs.Q and (0, 
so that Wy(L, t) ~ t K for times t <C L z , and saturates, 



Wr(L,t) ~ L a , for long times t L z , where a and 
z = &/k are the roughness and dynamic exponents of the 
interface defined by the local slopes, respectively. The growth 
equation for T(x, t), Eq.(|5}, may yield an interface T(x, t) 
that is not rough (a < 0). However, depending on the model 
symmetries, i.e., the form of the terms entering Q in (0, the 
corresponding local slope T(x, t) may turn out to be rough 
(q > 0), which straightforwardly implies anomalous scaling 
behavior of the surface h(x, t). 

We now analyze the dynamic scaling behavior of the sur- 
face slopes described by Eq.Q- Let us consider the expan- 
sion of Q(T) in powers of T, its derivatives, and combinations 
thereof (i.e. the leading-order gradient expansion) consistent 
with all the symmetries of the problem under study 12511 . The 
corresponding expansion may have several terms. However, 
the long wavelength limit (after renormalization) is governed 
by the most relevant term. We do not need to know explicitly 
the form of the most relevant term in the expansion, but we 
just assume that it generically behaves as Q(T) — > b~ n T m , 
with n > and m > 1, when the scale transformation 
x — > &x is applied |25|. Note that this term does not have 
necessarily to be present in that precise form in the expansion, 
but it can be generated by renormalization of other terms. To 
illustrate this, consider for instance the V(Vh) 3 nonlinearity 
that is known to renormalize to the linear V 2 h term 126ll . In 
other words, we may not know a priori the values of n and to, 
but we can certainly assume that such a scaling must exist. 

The dynamics of the surface slopes, Eq.Q, is always con- 
served by construction. This immediately leads to the hyper- 
scaling relation z = 2a+d+2 for the slope critical exponents, 
which is an exact outcome from the non-renormalization of 
the noise intensity D for systems with conserved dynamics 
lEM - This hyperscaling relation is exact, independent 
of the detailed form of Q, and allows us to directly link the 
anomalous time exponent to the dynamics of the local slopes 
by 2k = 1 - [d + 2)/z. 

Since we assume that the leading term of the expansion 
renormalizes as G(T) ~ b~ n T m , critical behavior can 
be obtained from dimensional analysis I23I EH Eill . Ap- 
plying the self-affine scale transformation x — > fox, t — » 
b z t, and T -» b & T to Eq.©, we obtain that the rele- 
vant terms scale as {31/ dt) -> b"- ? - (dT/dt), VG{T) -» 
^am-(n+i) \/g(^c^ a nd the conserved noise scales simply as 
T] c — ► b~ 1 ~ l - d+z ^ 2 rj c , since it does not renormalize. These 
three terms are equally relevant after rescaling only if the fol- 
lowing scaling relationships are satisfied, 

( (m - l)a + z - {n + 1) = 

\2a-z + d + 2 = 0, {) 

where we again obtain the hyperscaling relation due to the 
conserved dynamics of the surface slopes, as expected. These 
two basic scaling relations must be satisfied by the critical 
exponents of the surface slopes for any growth model with 
the appropriate (renormalized) values of n and to. From 
Q we can obtain formal expressions for the critical expo- 
nents of the local slope, a — (n — d — l)/(m + 1) and 
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K = a/z = (n - d - l)/[2n - d + (d + 2)m]. Although 
the values of n and to can only be determined by dynamical 
renormalization group techniques for each particular model, 
we do not need in this Letter to know their explicit values. 
For the sake of our argument, we only need to observe that 
there must exist constraints on the possible values that n and 
to can have if the local slopes T(x, i) have to be rough. The 
conditions d > and k > yield 

n>d+l, (8) 

which severely restrict the form in which the most relevant 
term of Q(T) can renormalize in surface models exhibiting 
anomalous roughening. Consider for instance 1 + 1 anoma- 
lous surface growth, the condition (|8j implies that the most 
relevant contribution in the expansion of the relaxation term 
must behave under rescaling as Q(Vh) — > b~ 3 (V h) m , where 
to can be any m > 1. The most relevant term of this type cor- 
responds to to = 1, e.g. the linear term V 4 h. In fact, the linear 
equation dth = — V 4 h + -q is known to exhibit (super-rough) 
anomalous scaling. However, even if the V 4 /i is not present 
(or not relevant), a most relevant nonlinear term might exist 
and satisfy Eq.(|8j, yielding a > and k > 0, and therefore 
anomalous scaling. 

In the following sections we analyze separately conserved 
and non-conserved surface dynamics and show that the re- 
quirement of local dynamics together with invariance of the 
growth equations under some symmetry transformations re- 
sults in the impossibility of intrinsic anomalous scaling. Cer- 
tain types of conserved dynamics, however, may exhibit 
super-roughening scaling properties. 

Conserved dynamics.- If the surface height h(x, t) is lo- 
cally conserved the growth equation © must obey the conti- 
nuity equation 

a t /7 = -V-j + 77(x,<), (9) 

where Q(Vh) = -V • j. The term j(Vh) = j eq + j neq 
includes equilibrium and nonequilibrium contributions to the 
current. Invariance under the choice of the height origin im- 
plies that the current j cannot depend explicitly on h 12511 . This 
type of equation has attracted much attention in recent years 
because of its relevance in the continuum theory of molec- 
ular beam epitaxial growth. We can follow the same pro- 
cedure as above to evaluate the scaling behavior of the sur- 
face h(x,i) from Eq.@ when Q(Vh) = — V • j. Conserva- 
tion implies that the noise scales as 77 — > 6~( d + z )/ 2 r\. For 
the velocity term we have (dh/dt) — > b a ~ z (dh/dt), and 
the most relevant term of deterministic part renormalizes as 
0(V7i) -> b am -^ n+m ^g{S/h). In order to have scaling the 
three terms must be equally relevant and we find 

f (m - l)a + z - (to + n) = 

\2a-z + d = 0. ( j 

These scaling relations are satisfied for any conserved model 
with the appropriate (renormalized) values of n and m. Re- 
garding the values of to and 77 the same comments as above 
apply here. 



We can eliminate z from Ea. dlQt to obtain (777 + 1)0; = 
777 + 77 — d. Equivalently, we have (m + l)d = 77 — 1 — d 
from Eq.0. This allows us to eliminate any dependence on 
the unknown exponents 77 and to and find the relationships 

a = d + 1 and z = z, (11) 

which are exact and valid for any growth model with con- 
served dynamics. The relations (II li immediately imply that 
for surface growth with conserved dynamics either the scaling 
is Family- Vicsek (d < and a < 1), or it is super-rough (if 
d > 0, then a > 1). In the latter case, we can calculate the 
local roughness exponent a; oc by making use of the scaling 
relation a\ oc = a — zk (valid whenever n, a > 0) to obtain 
that ai oc = a — a = 1, which shows that, if anomalous scal- 
ing occurs in conserved growth, it can only be in the form of 
super-roughening . 

Non-conserved dynamics.- In this case the lack of con- 
servation leads to the renormalization of the noise amplitude 
D and the violation of the hyperscaling relation in (II 01 . A 
one-loop dynamic renormalization group calculation for the 
spectral function (h(— k, t)h(k, t)) yields the flow of the noise 
intensity dD(l)/dl = D [z - d - 2a + g(l) 2 )], where g(l) 
is the coupling constant, which renormalization flow depends 
on the actual form of the most relevant nonlinearity. This 
shows that in the scaling relations (II 01 . only valid for con- 
served growth, we must replace the hyperscaling relation by 
z — d — 2a + g% = 0, where g» is the fixed point. This leads 
to the relation d = a — 1 + gt/{m + 1) instead of ( II It . which 
clearly shows that the naive scaling V/i ~ b 01 ^ 1 is generally 
incorrect for non-conserved dynamics. 

The scaling relations for the critical exponents of the sur- 
face slopes Q, and in particular the condition for anoma- 
lous scaling (|8}, are also valid in the case of non-conserved 
growth. Therefore, we can analyze the potential form of 
the nonlinear terms in Eq.@ that are required for anoma- 
lous roughening in the case of non-conserved growth. In 
general, the most relevant nonlinearity will be of the form 
Q(Vh) ~ (V" 1 ^) • • • (\7" N h). However, not any term can 
actually appear in a leading order expansion of Q. Under the 
assumption of relaxational dynamics 113 ill , the shift invariance 
symmetry h — > h + c greatly bounds the form of the possible 
nonlinearities. This was studied in detail by Hentschel in Ref. 
130], where he showed that for growth models obeying non- 
conserved relaxational dynamics the expansion of Q to lowest 
order compatible with shift invariance is generically given by 

d t h = 2T F'V 2 h + 2Y F"V[(Vh) 2 } ■ Vh - 

-sT[F -2F'{Vh) 2 ] +77, (12) 

where T and F are functions of (V/7) 2 , and primes denote 
derivatives with respect to the argument [30]. Ulterior expan- 
sion of T = T +T 1 (Vh) 2 + - ■ ■ andF = F Q +F 1 (Vh) 2 + - ■ ■ 
in powers of (V/7) 2 provides all the terms compatible with 
shift invariance. In general, all the possible nonlinear terms 
in the expansion (I12> have the form T 2a , (V • T)T 26 or 
T 2c VT 2 • T. This means that, at any order in the expansion 
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the only allowed terms compatible with shift invariance lead 
to n < 1, and these values are incompatible with the require- 
ment Eq. in any dimension d. Since any term compatible 
with shift invariance prevents condition (|8} from being satis- 
fied, we conclude that anomalous roughening (either intrinsic 
or super-rough) cannot take place in non-conserved growth 
equations. 

Conclusions.- We have shown that symmetries and con- 
servation severely bound the form of the most relevant terms 
that can exist in the long wavelength description of growth 
models exhibiting anomalous scaling. This has been exploited 
to show that intrinsic anomalous roughening cannot occur in 
local growth models, but some conserved dynamics may dis- 
play super-roughening for some types of interaction terms. 
Our conclusions are also valid in the case of growth driven 
by conserved noise in Eq.(0} after transforming d — » d + 2 in 
the scaling relationships (7) and (10), which leads to the same 
general conclusions. These results are useful in the search 
for continuum models and the general problem of assigning 
universality classes to experiments or discrete growth mod- 
els. Our results imply that disorder and/or non-local effects 
(like shadowing or bulk diffusion) must be responsible for the 
occurrence of intrinsic anomalous roughening in experimen- 
tal systems 

M OH El HI El. Regarding simulations 

of discrete models of local growth where anomalous scaling 
has been reported |6, 8], like for instance in simulations of the 
Wolf- Villain model 1 29] for ideal molecular-beam epitaxy, our 
results imply that those anomalous scaling exponents have to 
be taken with caution since they are effective exponents cor- 
responding to a transient, hence nonuniversal, regime. This 
conclusion is in agreement with most recent simulation re- 
sults of the Wolf- Villain model in large system sizes and long 
times iH and theoretical arguments l32il . 
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